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We discuss the cooling efhciency of ultracold Fermi-Bose mixtures in species-selective traps using 
a thermodynamical approach. The dynamics of evaporative cooling trajectories is analyzed in the 
specific case of bichromatic optical dipole traps also taking into account the effect of partial spatial 
overlap between the Fermi gas and the thermal component of the Bose gas. We show that large 
trapping frequency ratios between the Fermi and the Bose species allow for the achievement of a 
deeper Fermi degeneracy, consolidating in a thermodynamic setting earlier arguments based on more 
restrictive assumptions. In paxticular, we confirm that the minimum temperature of the mixture 
is obtained at the crossover between boson and fermion heat capacities, and that below such a 
temperature sympathetic cooling vanishes. When the effect of partial overlap is taken into account, 
optimal sympathetic cooling of the Fermi species may be achieved by properly tuning the relative 
trapping strength of the two species in a time-dependent fashion. Alternatively, the dimensionality 
of the trap in the final stage of cooling can be changed by increasing the confinement strength, 
which also results in a crossover of the heat capacities at deeper Fermi degeneracies. This technique 
may be extended to Fermi-Bose degenerate mixtures in optical lattices. 

PACS numbers: 03.75.Ss, 51.30.+i, 32.80.Pj, 67.60.-g 



I. INTRODUCTION 

Atomic physics and condensed matter physics now en- 
joy strong connections through the study of quantum 
transport in ultracold dilute gases Long-standing 
problems of condensed matter physics may be addressed 
by preparing controllable environments for the dynamics 
of cold atoms and by continuously tuning their interac- 
tions. This in turn allows for the study of fundamental 
features of high-temperature superconductivity using ul- 
tracold gases as controllable, analog computers of various 
model hamiltonians 

Degenerate Fermi gases were first produced in 1999 
Q, and more recently Fermi supcrfluid behaviour has 
been conclusively evidenced through the generation of 
vortices Q and the onset of critical velocities |^ in 
degenerate samples of ^Li. Weakly interacting Fermi 
gases arc difficult to bring to quantum degeneracy mainly 
due to fundamental obstacles in adapting cooling tech- 
niques successfully used for bosonic species. In par- 
ticular, the Pauli principle inhibits efficient evaporative 
cooling among identical fermions as they reach degener- 
acy. This issue has been circumvented by developing two 
cooling techniques, namely mutual evaporative cooling 
of fermions prepared in two different states and sympa- 
thetic cooling with a Bose species. In the case of dual 
evaporative cooling, a selective removal of the most ener- 
getic fermions in both the hyperfine states is performed. 
Provided that the initial number of atoms in each state is 
roughly the same, efficient dual evaporative cooling can 
be performed throughout the entire process. Limits to 



the minimum reachable absolute temperature using dual 
evaporative cooling have been addressed in resulting 
in a minimum reachable temperature T ~ ^i/ks, with 
/i the chemical potential of the Fermi gas (see also Q 
for a complementary analysis). Moreover, the number 
of available atoms Nf progressively decreases over time 
with a corresponding drop in the Fermi temperature Tp 
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proportional to Nj . The resulting gain in terms of 
a lower T /Tp degeneracy ratio is therefore limited, and 
the smaller clouds obtained at the end of the evaporative 
cooling are detrimental to detailed experimental inves- 
tigations requiring a large number of atoms, such as a 
quantitative mapping of the supcrfluid phases. In the 
case of sympathetic cooling using a Bose gas, the num- 
ber of fermions is instead kept constant, leaving aside 
losses due to background pressure and two- and three- 
body collisions, and the cooling efficiency depends on 
the optimization of Fermi and Bose collisional proper- 
ties, heat capacities, and, in the case of inhomogeneous 
samples, their spatial overlap. 

To date, the smallest Fermi degeneracy achieved with 
both cooling techniques is in the T/Tp > 5 x 10~^ range 
d, H]. This limitation has not precluded the study of 
temperature-independent features of degenerate Fermi 
gases, such as quantum phase transitions related to un- 
balanced spin populations [l^, El, [13, [111 or the effect 
of Fermi impurities in the coherence properties of a Bose 
gas [13, [iBl- However, the study of more conventional 
phase transitions in which the temperature is the key 
parameter is still uncharted territory and, as discussed 
for instance in [l^, [13, [3 1 this requires the achievement 
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of degeneracy factors T/Tp ~ 10~^ or lower. Unconven- 
tional pairing mechanisms that are unstable at higher 
T/T-p could then be observed, and the phase diagram of 
Fermi atoms in the degenerate regime could be mapped 
in a wider range of parameter space. Moreover, the study 
of ultracold Fermi-Bose mixtures is an interesting subject 
in itself, acting as the counterpart of the ■^He-^'He liquid 
mixtures extensively investigated at much higher densi- 
ties and temperatures. 

Considering the novel physical insights that deeper 
Fermi degenerate gases and Fermi-Bose mixtures may 
provide, it is relevant to discuss the limitations to reach- 
ing the lowest T/Tp in realistic settings available by 
means of sympathetic cooling, and ways to overcome 
them. Here, wc discuss two different techniques to over- 
come the apparent T/Tp ~ 10~^ limit observed so far, 
based on optimized heat capacity matching with species- 
selective traps or with lower dimensionality traps. The 
paper is organized as follows. In Section II, after briefly 
reviewing previous results on sympathetic cooling in 
spccies-sclcctive traps, we determine the time evolution 
of the temperature of the mixture in a particular class of 
species-selective traps through a thermodynamical anal- 
ysis, and we subsequently include the effect of the spa- 
tial overlap between the thermal component of the Bose 
gas and the Fermi gas. The main novelty with respect to 
previous semiquantitative analyses is that including both 
temporal and spatial dependence in the thermodynamics 
of sympathetic cooling leads to an optimization of the 
heat capacity provided that a time-dependent trapping 
frequency ratio is implemented. As an alternative to this 
optimization procedure, in Section III we discuss heat 
capacity matching resulting from lower effective dimen- 
sions for trapping, exploiting the strong dependence on 
dimensionality of the density of states of the Bose gas. A 
simpler protocol for optimizing heat capacity is available 
with nearly one-dimensional Fermi-Bose mixtures, which 
requires time modulation of the trapping strengths of the 
Fermi and Bose gases in the last cooling stage. Broader 
considerations on generic trapping settings are then dis- 
cussed in the conclusions. 



II. HEAT CAPACITY MATCHING THROUGH 
SPECIES-SELECTIVE TRAPPING 



After discussing the evidence for a correlation between 
the degeneracy factor T/Tp and the trapping frequency 
ratio between the Fermi and the Bose species, we review 
previous results on the use of species-selective traps and 
their limiting assumptions. We then relax these assump- 
tions with a thermodynamical analysis also including the 
effect of the partial overlap between the Fermi gas and 
the thermal component of the Bose gas in the specific 
case of bichromatic traps. 



Qualitative considerations on heat capacity 
matching 



Evaporative cooling has been instrumental in reaching 
Bose degeneracy for dilute atomic gases. Extensive anal- 
ysis have already addressed the dynamics of evaporative 
cooling of a Bose gas [l^ [13, [2TI. [2^ , using Monte-Carlo 
[H, US , mean- field analysis |25j | , and beyond [26t , in- 
cluding also more detailed effects [13, [H, H^. These 
studies have been also extended to the case of separate 
Bose and Fermi clouds (30j or Fermi-Bose mixtures using 
the quantum Boltzmann equation [sij or other semiclas- 
sical models [13 • Some generic features of the effective- 
ness in cooling fermions through a Bose gas can be ad- 
dressed based on the insights first discussed in [sl] and 
then analyzed in more detail in [sj, H^. The heat ca- 
pacity of a degenerate Fermi gas depends linearly on its 
temperature, being for a harmonically trapped gas equal 
to Cf ~ ir'^ksNfT /Tp, while a harmonically trapped 
degenerate Bose gas has a cubic dependence on temper- 
ature Cf, ~ 10.8 kBNb{T/Tc)^. The degeneracy param- 
eter can be written in terms of the ratio of heat capaci- 
ties: T/Tp. ~ 0.35(w6/w/)3/2(C6/C/)i/2 ^ By assum- 
ing that sympathetic cooling loses efficiency when the 
heat capacity of the Bose gas matches exactly that of 
the Fermi gas {Cb = Cf) wc obtain a conservative limit 
on the attainable T/Tp vs. ojf/uob space, depicted in 
Fig. 1 by the upper line. In the hypothesis that some 
residual cooling occurs when Cf, < C/, for instance with 
cooling stopping when Cf,/C/ < 0.1, we obtain the lower 
line in Fig. 1. Realistically, we do expect that sympa- 
thetic cooHng will be quenched when 0.1 < Cf,/C/ ~ 1, 
i.e. in the region delimited by the two lines. In Fig. 
1 we also plot the minimum T/Tp as obtained by the 
various running experiments with Fermi-Bose mixtures. 
Although the number of explored Fermi-Bose mixtures 
is limited and the diverse technical solutions for trap- 
ping and cooling may provide alternative explanations 
[36| . a correlation between the trapping frequency ratio 
ujf/ujb and the minimum achieved T/Tp seems corrob- 
orated by the actual results and invites more quantita- 
tive attention. Further analysis provided insight into the 
cooling limitations in different trapping conditions and 
for different species combinations [4l|. It was found that 
significant gains in T/Tp could be achieved for stronger 
relative Fermi-Bose confinements than the natural oj / /wf, 
provided by the mass ratio between the two species. In 
[4lj the focus was on an equilibrium situation at nearly 
zero temperature, and to develop a more comprehensive 
understanding of the cooling process a dynamical frame- 
work is needed in which crucial finite temperature ef- 
fects for the Bose gas are taken into account. Here we 
specialize the analysis of the thermodynamics as ^Li is 
sympathetically cooled by evaporating *^Rb in an optical 
dipole tra p. T his mixture, expected to optimize cooling 
42|, is currently used in various laborato- 
45[ both for studies of Fermi superfluidity 
and for the formation of ultracold molecules [4^, taking 
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FIG. 1: (Color online) Plot of experimentally obtained T /Tf 
versus the trapping frequency ratio ujf/iut- Lines indicate 
the theoretically predicted range of T/Tp values based on 
a semiquantitative heat capacity matching argument as dis- 



cussed in the text. The experimental data are taken from 
(Amsterdani)^^[i] (MIT)j3l (Florence), ^ (Boulder), 
(Houston), [sil (Paris), g3(Tubingen) . 



advantage of the large electric dipole moment of Li-Rb 
molecules (irj . 

In all previous discussions of species-selective traps a 
static picture was assumed for the cooling of the Fcrmi- 
Bose mixture, relying on the analysis already available in 
the case of single-species trapping [i^ . A critical element 
of this cooling is to balance the thermalization and loss 
rates; rapid trap modification allows for minimization of 
atom losses but is limited by the requirement that wc 
go proceed slowly enough to keep the system at thermal 
equilibrium. It has been shown [49[ that thermodynamic 
equilibrium is maintained if the ratio between the poten- 
tial depth and the atomic cloud temperature, given by 
r] = AU/kBT, is kept constant. The time-dependence of 
the potential depth is then given by (isj : 



AUjt) 
AC/, 



t 

1 + - 
r 



(1) 



where AUi is the initial potential depth, r]e = —2(77' — 
3)/r]', and t^^ = [2/i)ri'{T] - 4) exp(-77)7,, with -q' = 
77 + (j? — 5)7(77 — 4) and 7,; the initial elastic collision 
rate. With the time dependence of A[/ determined in 
this way wc can obtain the other relevant quantities in 
the process (number of particles, temperature, density, 
scattering rates) resulting in scaling laws similar to Eq. 
([1]). A potential depth/ temperature ratio of 7;=5 to 10 
is considered to yield optimal efficiency in evaporative 
cooling [4^ . This approach can be applied to either the 
simple case of evaporative cooling of a single Bose species, 
or the sympathetic cooling in a Fermi-Bose mixture. In 
both cases, Eq. ([1]) describes the potential depth of the 
Bose species AC/f,, with the implicit assumption that the 
presence of the fermions does not drastically affect the 



evaporation and cooling of the bosons. This is justified 
since all dual-species systems trap at least an order of 
magnitude more bosons than fermions, but this approxi- 
mation may suffer towards the end of the cooling process 
when a majority of the bosons have been evaporated. 

In the case of a bichromatic optical dipole trap, from 
AUb wc determine the required power Pi of the laser 
confining both species, and then for a targeted trapping 
strength ratio Wf /lOh we can calculate the required power 
P2 for the Bose-deconfining laser, and finally determine 
the fermion potential depth AUf. In this way we have in- 
dependent control of the spatial size and potential depth 
for each of the two species, allowing us to either main- 
tain a constant ujf/ujb throughout cooling or adjust the 
relative trapping strengths during the process. 

The exact way in which the temperature is determined 
from the trapping parameters depends upon both the 
model and the proposed cooling strategy. Eq. ([1]) merely 
identifies a limit to the cooling rate if thermodynamic 
equilibrium is to be maintained throughout the process, 
and assuming a constant ratio 77, as it was discussed in 
[3^ . In practice, these assumptions can be relaxed by 
using a more dynamical model based on energy balance, 
as we discuss in the next section. 



B. Thermodynamical balance 

To move beyond the limiting assumptions in [33 |. we 
start our analysis by considering a dual-species system 
at thermodynamic equilibrium. The trap parameters are 
then suddenly changed in order to force some bosons 
to evaporate; then one waits for a new thermodynamic 
equilibrium before applying another evaporation step, in 
analogy to the scheme discussed in [s^l ■ The step-by-step 
temperature reached in this way is determined by energy 
conservation for the Fermi and Bose gases [5l| . 

For concrctcness, we consider a species-selective trap- 
ping scheme as described in (s^ . A mixture of Nf 
fermions and Nh bosons is confined into a bichromatic 
optical dipole trap tailored by two lasers of wavelengths 
Ai, A2 and powers Pi, P2- In order to obtain quasi an- 
alytic results, we will approximate the trap potential by 
a truncated harmonic potential properly reproducing the 
bottom curvature and the depth of the well. This ap- 
proximation becomes exact for energies small with re- 
spect to the trap depth, a condition satisfied in the cases 
discussed below. The presence of the second laser allows 
one to make the trap parameters of the fermion species, 
namely the characteristic frequencies iUfx,ujfy,u!fz and 
the depth AUj, different from the corresponding bo- 
son parameters, ujbx,^by,^bz and AUb- As a conse- 
quence, the ratio ujf/ujb, for each species s = 6, / we 
define lUs = {tUsx^^syi-^szy^^ , can be varied from its mass- 
determined value for P2/P1 = to an ideally arbitrary 
large value for P2/P1 approaching a positive critical value 
[54I . The evaporation steps are carried out by decreasing 
the power Pi of the reference laser while maintaining the 
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ratio P2/P1 at a constant value. In this way, the trap 
depths AUb, AU-f, which are proportional to Pi, and the 
frequencies w;,, (proportional to y/Pi), decrease while 
the ratio ujf/ujb remains constant [53[. 

At the end of an evaporation step in which the power 
of the reference laser is changed from Pi to Pi + dPi 
and once thermodynamic equilibrium is reestablished, 
the temperature changes from T to T + dT according 
to an energy balance equation of the form 



{AUb + ^keT) dN^"" = dEb + dE 



(2) 



where dN^^ is the number of bosons in the excited states 
at temperature T that evaporate and AUb + ^ksT (with 
< ^ < 1) is the mean energy per evaporated boson [2^. 
The quantities dEb and dEf are the energy changes of the 
trapped boson and fermion species, both at temperature 
T. 

By observing that due to a change of Pi , all the quan- 



tities T, LOb, AUb, ^f, AUf change and using for iV^^, Eb 
and Ef the expressions provided by Eqs. (|A9p . (|A8[) and 
(|A12p . we have 



dEb 
dEf 



dT ^ 



f -TT^diUb 
OUJb 

dEb, 
t; — dujb - 



(5) 



It is worth to point out that, within the approximation 
used, Ef does not depend on AUf and the number of 
fermions Nf remains constant. By evaluating the above 
partial derivatives and inserting the result into Eq. 
we arrive at 




where we used diT/Tp) = dT/Tp - [T/TF)dTF/TF and 
dTp/Tp = dPi/{2Pi), which stems from the proportion- 
ality of Tp to Ljf, namely ksTp — {(SN fY/'^hujf, see Eq. 

(|XT3l) . Like iOf/uJb, the value of AUbPi^''^ /kBTp is a 
constant determined by the value of the ratio PijPi, and 
p(a;) and q(x) are defined as 

p{x) = / ^^dt + i^-, (7) 
Jo - 1 - 1 

q(x) = 3x / dt-4: / dt 

Jo e* - 1 7o e* - 1 

Note that p(0) — q{Q) = whereas p{x) ~ 6^(4) and 
q{x) ~ 6C(3)a: - 24C(4) + 6C(3)^ for x > 1, C being the 
Riemann zeta function. 

Equation ([S]) is a nonlinear ordinary differential equa- 
tion which allows for the determination of T/Tp{Pi) dur- 
ing the evaporative cooling. Observing that both func- 
tions p{x) and q{x) are non negative for a; > and as- 
suming for simplicity ^ = 0, the qualitative behavior of 
T/Tf{Pi) is as follows. The numerator of the last frac- 
tion in Eq. ^ is always positive. If we start from ini- 
tial values of Pi and T/Tp respectively not too smaU 
and not too large, the argument of the functions p and 
q is large with respect to unity, which is equivalent to 
state that AUb ^ ksT, a fact that also justifies the 



choice ^ = 0. The denominator of the last fraction in 
Eq. ^ is thus also positive so that T jTp decreases by 

decreasing Pi . The decrease may be faster or slower than 

1/2 

Pi depending on the value of the constants lo f /tob and 

— 1/2 

AUbPi /kBTp. Eventually, however, the last fraction 
in Eq. ([5]) becomes larger than 1 so that a second regime 
starts in which T/Tp decreases faster and faster. As a 
consequence, the argument of the functions p and q de- 
creases and the denominator of the last fraction in Eq. 
([6| approaches 0. A singular point is thus reached in 
which d{T/Tp)/dPi = 00 and T/Tp > 0. A numeri- 
cal study also shows that Eq. ^ has a discontinuity 
at the singular point with unphysical negative values of 
T/Tp on the left. The value of T/Tp at the right of the 
singular point represents the minimum achievable T/Tp 
during the cooling, provided that the initial number of 
bosons is sufficiently large so that they are not completely 
evaporated before the singular point is reached. 

The behavior of T/Tp as a function of Pi is shown in 
Fig. 2 for different initial conditions and different values 
of LOf/ojb- We stress that whereas the minimum T/Tp de- 
pends very little on the details of the cooling, the number 
of atoms of both species and initial values of tempera- 
ture and reference laser power, we observe a substantial 
decrease of the minimum achievable T/Tp, by increasing 
the ratio lo / /ub- This is in agreement with a previous pre- 
diction based on a rough matching of boson and fermion 
heat capacities [sJl. In fact, the singular point of Eq. ^ 
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FIG. 2: (Color online) Dependence of the degeneracy factor 
T/Tf upon the confining laser power Pi during sympathetic 
forced evaporative cooling as determined by Eq. (|6]). The 
system is a mixture with Nf atoms of ^Li and Ni, atoms of 
^^Rh trapped in a bichromatic optical dipole trap shaped by 
two lasers of power Pi and P2 at the wavelengths of Ai = 
1064 nm and A2 = 740 nm for the ^Li-^'^Rb mixture as cho- 
sen in [3^. Two sets of curves are shown for different initial 
conditions and for different values of the ratio P2/P1, kept 
constant during the evaporation. The fermion-to-boson trap- 
ping frequency ratio, determined by P2/P1, is ujf/uib ~ 2.443 
for P2/P1 = 0, ujf/ujb = 8.186 for P2/P1 = 0.18 and 
LVf/uJi = 15.911 for P2/P1 = 0.23. For the same P2/P1 values, 
the other constant AUbPi ^'''^ /ksTp which appears in Eq. (j6]) 
amounts to 18.73, 3.70 and 0.59, respectively. For simplic- 
ity, we set ^ — 0. We assume Nf = is constant during 
the evaporation and, for the continuous (dashed) curves the 
initial number of bosons is Nb ^ 2 X 10'' {Nb = 2 X 10^), 
Fermi degeneracy T/Tf = 1 {T/Tf = 10) and Pi = 1 W 
(Pi — 5 W). The minimum achievable T/Tp, corresponding 
mathematically to a singularity of Eq. ([Sjl and physically to 
a fermion-boson heat capacity equality, (a) does not depend 
on Nf and A^i,, provided that A^;, is sufficiently large and (b) 
depends only slightly on the initial conditions for T/Tf and 
Pi, but (c) decreases appreciably if the trapping frequency 
ratio LOf/ujb is increased. 



is defined by the condition 

fATT ^ci. i^^^^b' dE, dEf 



(9) 



For AUb ^ ksT, a condition well satisfied at the singular 
point, we have 



1 dEb 6C(3) 
dT ksT dT 24C(4) 

therefore Eq. ([9]) is equivalent to 

AUb 



0.28 



kBT 



(10) 



(11) 



FIG. 3: (Color online) Dependence of the degeneracy factor 
T/Tf upon the number of bosons (normalized to the number 
of fermions) during a sympathetic forced evaporative cool- 
ing driven by the laser power Pi as determined by Eq. 
The parameters are the same as in the case of the continu- 
ous curves in Fig. 2, with the initial values of A^i, = 2 x 10^, 
T/Tf = 1, and Pi = 1 W". It is evident that the use of 
larger Uf/ujb ratios allows us to reach a deeper Fermi degen- 
erate regime, which amounts to a gain by almost two orders 
of magnitude difference in the case of the larger uof/uob ratio. 



The sharp drop observed in T/Tp before the singu- 
lar point deserves some comments. At temperatures suf- 
ficiently low with respect to Tp and Tc, the energy of 
the mixture E^ -\- Ef is dominated by the zero temper- 
ature Fermi energy Ef{Nf,0), see Appendix A. In this 
case, the right-hand side of the energy balance equation 
([2|) can be approximated by {dEf{Nf,0)/dujf)dujf, and 
in the proximity of the singular point (i.e. at low tem- 
peratures), following Eq. (fTTj) the left-hand side of the 
same balance can be written as {'jCb — Cf)dT, where 
7 ~ 0.28AUb/kBT — 1. The approximated energy bal- 
ance thus gives 



dcjf iNfkeTp 
ujf jCb-Cf 



(12) 



where Cs = dE,/dT with s^bj. 



From Eq. (|T2|) we see that a small decrease of the fermion 
trapping frequency induces a temperature decrease, the 
size of which depends on the value of N jksTp / {-^Cb — 
Cf). Note that the denominator of this ratio contains a 
difference, not a sum, of the boson and fermion specific 
heats. The divergence of the derivative dT/dcuf predicted 
at the critical point is certainly unphysical: we expect 
that a long time is needed to re-equilibrate the system in 
a freezing step T ^ T + dT with dT large. In this case, 
dissipative phenomena should be taken into account by a 
more complicated model in which the singular point will 
be substituted by a minimum. However, this does not 
change the meaning of the lowest reachable T/Tp which 
is the point where fermion and boson specific heats do 
match. 

In Fig. 3 we plot the dependence of T/Tp upon the 
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number of bosons, normalized to the fermion number (as- 
sumed to be constant during the evaporation process). It 
is manifest that deeper Fermi degeneracy factors are ob- 
tained for higher trapping frequency ratios. This plot has 
to be compared to the one presented in Fig. 2 of in 
which T/Tp was shown versus a similar quantity (in our 

notation (A''^'"' — Nb)/Nf where iV^'"' is the number of 
initial bosons prior to evaporation). While our analysis 
confirms that the initial decrease in T/Tp is faster for 
lower LUf/uji, ratios, thus suggesting more cfhcicnt cool- 
ing - if measured by the drop in T/Tp per unit of boson 
removed in the evaporation process - we also notice that 
the evaporation process extends much further for larger 
LOf /cob and allows to reach deeper Fermi degeneracy fac- 
tors before stopping. Therefore, if the goal of the cooling 
process is the achievement of the lowest T/Tp degeneracy 
parameters rather than saving bosons during the evapo- 
rative process |54|. the use of larger uj f juj^ is beneficial, 
at variance with [32l |. 

The main issue in [3^ is that no temporal depen- 
dence for the trapping frequencies was assumed (unlike in 
Eqns. (HHSl)); which is unphysical for any realistic evap- 
orative cooling strategy involving species-selective trap- 
ping strengths. The discussion in [33| also struggles with 
issues arising from more practical limitations, as com- 
mented in [4l|, in particular the fact that all the Fermi- 
Bose species available in practice will be affected by issues 
of spatial overlap in usual confining potentials, due to the 
smaller mass of the Fermi species (apart from the never 
considered '**'K-^'^Na mixture). Also, for Wb/w/ > 1 the 
superfluid critical velocity of the Bose gas will become 
larger than the Fermi velocity, inhibiting scattering be- 
tween fermions and bosons and then sympathetic cool- 
ing [55I . [5^. Finally, a weaker trapping frequency for 
the Fermi species corresponds to a lower potential en- 
ergy depth with respect to that of the bosons, resulting 
in significant fermion losses during the forced evaporative 
process of the bosons. As we discuss in the following, a 
large ujf/ujb ratio is also beneficial in terms of improv- 
ing the spatial overlap between fermions and bosons and 
therefore the cooling efficiency. 

C. Spatial overlap between the fermion and the 
thermal Bose atoms 

An even more realistic analysis of the cooling dynamics 
must also account for the intrinsic inhomogencous char- 
acter of the trapping potential as this will result in in- 
complete overlap between the Fermi and the Bose gases 
and consequently a decreased cooling efficiency. As intro- 
duced in |4l| , we express the spatial overlap between 
two clouds of densities pi and pj, as 

e,, = {N^N,)-'/^ J p\l\r)p'l\T)d\ , (13) 

where i,j refer to the fermion (F), Bose condensate (B) 
or thermal boson (T) density profiles. The fraction of 



atoms that share the same region of space is thus given 
by e^^ . The Bose atoms available for cooling are those 
having a non-zero specific heat, i.e. those in a thermal 
state. Finite overlap between the thermal bosons and the 
Fermi atoms will result in a decrease in the cooling rate 
with respect to the case of ideal overlap. An accurate 
evaluation of the cooling rate should take into account 
kinetic equations for the two interacting gases. In a pes- 
simistic, conservative fashion, we can assume that the 
cooling rate q is decreased by a factor equal to the frac- 
tion of atoms that can actually exchange energy without 
any mass transport involved, as (jcooi — ^ Cft^cooI- The 
minimum attainable degeneracy parameter T/Tp corre- 
spondingly increases as T/Tp T/Tp. This static 
estimate does not take into account the timescale over 
which fermions and bosons exchange energy through clas- 
tic collisions, and the particle relocation along the trap 
volume, but it can be considered as an upper limit to 
the effect of partial overlap. This analysis requires the 
density profiles of the condensate and the non-condensed 
thermal boson to computed independently, following the 
discussion of a Fermi-Bose mixture at finite temperature 
reported in (57| . 

The dependence of the spatial overlap parameters Efb 
and Eft on temperature, boson number, and trapping 
frequency ratio is shown in Fig. 4. As the temperature 
drops below Tc a finite condensate fraction appears; the 
fermion-thermal boson overlap starts to decrease, and the 
fermion-condensate overlap increases. There is a limit to 
Efb for a given ujf/ujb, however, since the cloud of less 
massive fermions will have a much larger spatial radius 
and thus a strong relative confinement is required to im- 
prove the fermion-Bose condensate overlap. The overlap 
dependence on number ratio is rather straightforward, 
with a gradual decrease in fermion-condensate overlap as 
bosons are evaporated and an almost flat behaviour for 
Eft since N^'^ is roughly constant for a given temperature 
and boson losses manifest as a decrease in the condensate 
fraction and have minimal effect on the thermal cloud. 

The effect of the spatial overlap on the minimum reach- 
able Fermi degeneracy factor T/Tp is depicted in Fig. 5, 
with the choice of interspecies scattering length of the 
Li-Rb system corresponding to the pessimistic scenario 
of repulsive interaction (see for the issues related 
to its measurement). It is evident that about one or- 
der of magnitude may be lost in the achievable minimum 
degeneracy factor when the overlap factor is taken into 
account, although we conjecture that with a full kinetic 
analysis the actual result will be located in between the 
two curves. The situation is analyzed in more detail in 
Fig. 6, which shows Eft versus cuf /cub and T/Tc- During 
the cooling process, i.e. as T/Tc decreases, the opti- 
mal overlap is shifted to larger values of ujf /tUb, until the 
minimum of the Fermi degeneracy is reached as shown 
in the contour plot of Fig. 7. This suggests the use of 
a time-variable trapping strategy, initial with lower val- 
ues of ujf/ujb, then increased in time by increasing the 
power ratio P2/P1. Such a time-dependent relative con- 
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FIG. 4; (Color online) Overlap parameters efb (dashed line) and epT (solid line) versus (a) T/Tc, (b) Nt/Nf, and (c) Uf/coi,. 
Plots made with other two parameters held at given values of Nt/Nf — 15, T/Tc— 0.3, ujf/ub=3, in the case of a scattering 
length for Rb of ai,j,=+5.8 nm and an interspecies scattering length a/i,=+0.5 nm (with the same values also used for the plots 
in Figs. 5-7). 




FIG. 5: (Color online) Plots of the minimum T /Tf reached 
during the evaporative cooling process in an optical dipole 
trap, without (line) and with (dots) the fermion-thermal spa- 
tial overlap taken into account. The curve and dots start at 
the lowest possible value of Uf/uji,— 2.443, which is the 'natu- 
ral' trapping frequency ratio for a ®Li-®''Rb mixture with the 
deconfining laser beam switched off (P2=0). 



fincmcnt strategy is not the only way to optimize the 
Fermi degeneracy factor, however. In the following, we 
will discuss a similar procedure which exploits the advan- 
tages of reducing the dimensionality of the system when 
Fermi degeneracy is approached. 




a 4 S S 10 ]2 



FIG. 6: (Color online) Three-dimensional plot and two- 
dimensional contour plot of the Eft dependence on ujf/ui,, 
and T /Tc- Significant overlap values for trapping ratios in 
the range oi ujf /uj^ ~ 3-7 are evident at relatively large tem- 
perature ratios T /Tc. At lower temperatures, the optimal 
overlap is achieved at higher trapping frequencies ratios. The 
optimal path maximizing the overlap is highlighted by the 
dashed line in both plots. 



III. HEAT CAPACITY MATCHING THROUGH 
LOWER DIMENSIONALITY 

As an alternative to the cooling strategy described 
above, we discuss here the possibility to match the heat 
capacity of Bose and Fermi gases at the lowest possi- 
ble T/Tp by exploiting lower dimensionality traps. Ul- 
timately, the mismatching between the specific heats of 
Bose and Fermi degenerate gases depend on the scaling 
of the heat capacities with temperature, and this in turn 
depends upon the dimensionality of the Bose gas. As 



discussed in [sol . [60j and demonstrated in [6l|, [B^I , a dra- 
matic increase in the trapping frequency in one (or two) 
trapping axes will result in an effective two- (or one-) di- 
mensional system. This in turn allows for a better match- 
ing of the heat capacities since the Bose gas dependence 
on temperature will become milder than in the full 3D 
case. In order to gain quantitative insights on how to 
realize such a matching, we first consider noninteracting 
gases in a harmonic potential, with the mmrber iV/.t, of 
particles fixed. 
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FIG. 7: (Color online) Contour plot of the Fermi degeneracy 
factor T/Tf versus the Bose degeneracy factor T/Tc and the 
trapping frequency ration u / /uJb ■ The lower Fermi degeneracy 
factor of T/Tf — 0.02 is obtained for iJf/uJb — 8. 

The total number of fermions and bosons is evaluated 
as (+ for fermions. — for bosons): 

oo 

NfA^^^T) = Y. ,iE,J/lT^, (14) 

where gj is the degeneracy of energy level Ej , fi the chem- 
ical potential, ks Boltzmann constant, and T the tem- 
perature, with the number of particles Nf j, fixed. Solv- 
ing this equation numerically for fi = ^^{T), we can then 
calculate the total energy 

oo p 

and from this we obtain the heat capacity as C(T) — 
dE/dT. In the calculation below, we assign Nf^b = 10^, 
and assume the initial trap to be isotropic with the trap- 
ping frequency w = 27r x 15.87 kHz. The numerical calcu- 
lations for 2D and 3D are straightforward, while for ID 
some approximations are necessary to reduce the simu- 
lation time to realistic values. 



A. 2D and 3D traps 



FIG. 8: (Color online) Heat capacity curves of bosons and 
fermions in 3D (dashed curves) and 2D (solid curves) in an 
isotropic harmonic trap, with the bosons exhibiting nonmono- 
tonic behaviour. The crossing point between the Bose and 
Fermi curves for the 2D case occurs at a slightly higher T/Tf 
value, ruling out its use for a more favourable cooling of 
fermions. 

The upper limit Q in the summation should be infinity 
in principle, but a value of Q=1500 is sufficient for nu- 
merical convergence. For the given parameters -/V/.6 and 
Lo we solve the above two equations for ^ at different 
temperatures, and then calculate the heat capacities. As 
depicted in Fig. 8, the Bose and Fermi heat capacities 
intersect each other at T ~ 0.293rF, with the three- 
dimensional Fermi temperature Tp^ — {QN f Y/'^kuj /ks- 
The situation is very similar for atoms trapped in a 2D 
isotropic harmonic potential, except that now the en- 
ergy eigenvalues are given hyE'j^ = {j + \)fujj, with de- 
generacy Qj = j -\- \, and Q needs to be increased to 
~ 10^ to achieve adequate convergence. The two heat 
capacities curves intersect each other at T ~ 0.308 Tp 
where the 2D Fermi temperature Tp is given by T]P = 
{2N fY/'^hoj/kB- From Fig. 8 we see that going from a 
full 3D to a 2D system actually slightly worsens the heat 
capacity matching, yielding a higher T/Tp at the point 
where Cf, and C/ intersect each other. However, further 
reduction to a ID system results in complete matching 
of the heat capacities, as we will see below. 



For atoms trapped in a three dimensional harmonic 
potential V{r) = muf^r^ /2, the energy eigenvalues Ej 
(j=0,l,2,. . . ) are given by = {j + 3/2)fiw. Since the 
trap is three-dimensional and isotropic, the degeneracy 
gj of the energy levels is given by gj = (j + l)(j -I- 2)/2. 
The number of particles for bosons (— ) and fermions (+) 



N 



3D 



1 

2 ^ e 



(j- + l)(j + 2) 



{{j+3/2]nuj-p.)/kBT _|_ I 



(16) 



B. ID trap 

For atoms in a ID trap, the energy eigenvalues are 
given by Ej^ ~ (j + l/2)huj, with degeneracy gj = 1. 
If we follow the same steps as the 2D and 3D cases, the 
upper summation limit Q (see Eq. p6p ) must be ex- 
tremely large in order to reach convergence and approx- 
imations are required for the one-dimensional trap anal- 
ysis. These approximations are outlined in Appendix B 
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fermions by changing the ratio of the two trapping fre- 
quencies in a 2D trap. Thus one possible solution to 
improve the cooling efficiency is to first evaporate in a 
3D trap and then, when the Fermi degeneracy starts to 
reach about T/Tp — 0.3, to increase the trapping fre- 
quencies achieving a quasi-one dimensional system, then 
continuing the evaporation process. A possible limita- 
tion of this technique comes from the larger coUisional 
loss rate as a result of the increased confinement. Also, 
as studied for achieving Bose condensation of hydrogen 
atoms, the nearly ID character of the evaporative cooling 
[6^ may lead to non-ergodic evaporation limiting its effi- 
ciency [66[ , although this has not prevented achievement 
of Bose degeneracy [63] ■ 



FIG. 9: (Color online) Heat capacity curves of bosons (black, 
dashed) and fermions (red, continuous). We consider for 
simplicity the case of equal number of bosons and fermions 
N^, = Nf = 10* with equal mass mj, — nif. Curves in plots 
(a), (b), (c) and (d) are in 2D traps with k — coy/ujx equals to 
1, 2.5 X 10^, 2 X 10*, and 5 x 10* respectively, with the case 
in (d) showing an initial plateau of the heat capacities at ~ 1 
before reaching the high temperature limit of a 2 D gas. 

and with them we obtain the Fermi and Bose energies, 
from which the heat capacities are evaluated. To take 
full advantage of the gain in the heat capacity matching 
in ID traps it is helpful to investigate how heat capacities 
of bosons and fermions change when we gradually reduce 
the dimension of the trap. We assume that the atoms 
are first trapped in a 2D trap with trapping frequen- 
cies along the two dimensions tu and kto, and then this 
relative confinement parameter k is gradually increased 
from 1 to infinity. The system will be effectively ID when 
ksT <c khu. We evaluate the heat capacities of bosons 
and fermions at different values of k. The results are 
shown in Fig. 9 (a-d). As k increases, the shape of the 
heat capacity curve of bosons becomes more similar to 
that of fermions, in the sense that it slowly loses the 
peak structure and the curvature near zero temperature 
begins to resemble that of fermions. If k is further in- 
creased (case (c)) there will be a region where the two 
curves completely coincide with each other. This is con- 
sistent with the previous result we obtained for an ideal 
ID trap. In Fig. 9 (d) the heat capacity curves of bosons 
and fermions are shown at an even higher aspect ratio, 
to emphasize the crossover from the 1 D case to the 2 
D case at high temperatures. In the ideal ID case, the 
heat capacity curves are identical, as s imp ly explained 
in the canonical ensemble approach [6^, [gJ]. Indeed, the 
total internal energy of fermions in a ID harmonic trap 
only differs from that of the bosons by the Fermi zero- 
point energy Eq = Nf{Nf — l)fkv/2, and therefore the 
two systems have identical heat capacities. 

The existence of this crossover indicates that we 
can control the heat capacity matching of bosons and 



IV. CONCLUSIONS 

We have examined the thermodynamics of evapora- 
tive and sympathetic cooling in a Fcrmi-Bosc mixture, 
and identified possible ways to achieve a lower Fermi de- 
generacy factor T/Tp. Thermodynamical considerations 
are based on general assumptions and measurable, phe- 
nomenological inputs, like heating rate and specific heat, 
and provide a solid framework to discuss cooling dynam- 
ics regardless of sophisticated microscopic models [68j . 
They also allow for comparison with ex per imental results, 
such as those discussed in [H, [zO, [zH, lli, Izl, Izii or for 
the inclusion of more realistic inputs such as the specific 
heat of an interacting gas . Although we have focused 
the attention on the particular ^Li-^^Rb mixture, the ex- 
tension to other Fermi-Bose combinations is straightfor- 
ward, furthermore benefiting from more favourable in- 
terspecies thermalization properties with respect to this 
particular mixture, for which limitations in clastic scat- 
tering and use of Feshbach resonances have been experi- 
mentally evidenced [i^, [tI] . 

Among the main results we have obtained, we have 
discussed two different cooling strategies: constructing a 
spccics-sclcctivc trap with independently tunable Fermi 
and Bose trapping frequencies, and creating traps with 
reduced dimensionality in the latest stage of evapora- 
tion. We have shown that different trapping ratios lead 
to distinctly different cooling trajectories. However, in- 
complete spatial overlap will not only result in a longer 
cooling time needed to attain a given temperature, but 
will also increase the temperature at with the heating 
rate will balance the cooling rate. When the progres- 
sive depletion of the bosonic thermal cloud is taken into 
account, optimized cooling requires time-dependent trap- 
ping strengths. Additionally, we have discussed how to 
exploit the strong dependence of the bosonic specific heat 
upon dimensionality to create nearly one-dimensional 
traps in the ultimate stage of sympathetic cooling. This 
will be of particular relevance for various planned studies 
of Fermi gases in optical lattices fll\ in which bichromatic 
optical traps arc not viable. 
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APPENDIX A: BOSE AND FERMI GASES IN A 
3D HARMONIC TRAP OF FINITE DEPTH 

We report here the expressions for the energy and the 
number of particles of ideal degenerate Bose and Fermi 
gases confined a in a three-dimensional harmonic trap of 
finite depth. 

For the trapped bosons wc start with the exact expres- 
sions of the energy and of the number of particles in a 
general 3D harmonic potential: 



oo oo 

E E 

ria;— riy—O 7iz—0 



E 



Er, 



CXD OO OO 



^^-E E E 



1 



n^—O Uy—O n~— 



(Al) 



(A2) 



where £;„^,„^,„, 
hwbzin^ + 1/2). 



fi^bxin^ + 1/2) + huJbyiny + 1/2) + 



As usual, wc write Nt 



y 



where iV° = (e^o.o.o/'=B^— l)^i is the number of bosons in 
the ground state and the number of those thermally 
excited. In order that the number of particles remains 
positive, it is necessary for the chemical potential to sat- 
isfy fJ,{T) < Eofl.o- At temperatures T < Tc, where Tc is 
the Bose-Einstcin condensation critical temperature, the 
chemical potential is frozen to its maximum value. In 
general, fJ-{T) and all the other thermodynamic quanti- 
ties can be evaluated explicitly as a power series expan- 
sion in the two parameters fi{T) — £'o,o.o smd hajb/ksT, 
where ujb ~ {uJbx^jJbytjJbzY^^ [78|. In the case T < Tc and 
hujb <^ ksT, which is relevant to the experimental situa- 
tions discussed here, we can restrict to the lowest order 
and write 



Eb = 3C(4) 



= c(3) 



{nwbf 



hujh 



(A3) 



(A4) 



where C is the Riemann zeta function. Alternatively, the 
above two results can be obtained using the semiclassical 
density of states 



Qb{E) 



E^ 



dE hUbxflUJbyhuJbz "^ifl^bY 

and the continuum approximation 



(A5) 



Eb 



f 

Jo 



E 



gb{E)dE, (A6) 



Qb{E)dE. 



(A7) 



For a trap of finite depth schematized as a harmonic po- 
tential truncated at energy A[/f,, we thus write 



Eb = 



A (76 



E 



Qb{E)dE 



{ksTf 1 
(h^bf 2 



t3 



- 1 



(A8) 



AUb 



1 



E/ksT _ I 

[ksTf 1 t 
{hujbf 2 





3 



Qb{E)dE 



-dt. 



1 



(A9) 



Consider now a system of Nf fcrmions confined 
by a harmonic trap having characteristic frequencies 
ujfx,ujfy,u!fz- Under the condition hujfx,fiujfy,hujfz ^ 
ksT, a continuum approximation holds as in the boson 
case, so that the fermion counterparts of Eqs. (|Aip and 
()A2p can be simplified to 



Ef = 



Nf 



E 



4E-tj.)/ki 



^—gf{E)dE, (AlO) 



[ ei-'^)lsT^^ SfWE, (All) 



where Qf{E) = E'^ /2{hu!fY and w/ = {ufxUJfyUjfz)^^^ . 
Note that in the above expressions we have neglected 
the zero-point energy which is justified since for fermions 
fi 3> Eofi^o- The accuracy of the continuum approxima- 
tion introduced above is discussed in detail in [79]. The 
chemical potential fi(T) can be eliminated between Eqs. 
(|A10|) and (|A11[) by means of the standard Sommerfeld 
expansion in powers of temperature. By keeping terms 
only to second order in T, we have the well known result 

Ef{Nf,T) ^ EfiNf,0) + ^{kBTfgiEp), (A12) 
where the Fermi energy Ep is related to Nf by 



(A13) 



and the zero temperature term Ef(Nf, 0) is given by 

EfiNf,0)^j"'E,,iE)dE = ^. (A14) 

Making use of the relationship between the Fermi energy 
and the number of fermions, we finally obtain 

^ ^6'/^Ny\ujf + ^6-i/3iV^/^^%^. (A15) 
4 ■' 2 ■' 



Equations (jA13|) and (jAlSP are valid also for fermions 
trapped into a harmonic potential truncated at energy 
AC//, provided that AUj > Ep^ksT. 
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APPENDIX B: INTERNAL ENERGY FOR A ID 
TRAPPED IDEAL FERMI GAS 

The total atom numbers in one dimension N^^ and 
Nj^ are given by Eq. pi)) . with gj = 1. However, the 
upper summation hmit Q (see Eqs. (fT6)) ) needed for con- 
vergence becomes unreasonably large in one dimension 
jsO], and thus approximations are required. 

For fermions, we can directly replace the summation 
with an integral in the case of fc^T ^ fko: 



Having fixed N, we solve this equation numerically for 
fi at different temperatures to obtain fj, = n{T), and then 
apply the Euler-Maclaurin formula again to obtain the 
total energy E^^ . 

The discussion above holds for a one dimensional Bose 
gas. In the case of a trapping frequency range for which 
there is a smooth crossover from two dimensions to one 
dimension, we recall that the number of particles in a 2D 
trap is given by: 



E 

1 

huj 



1 



p{U + l/2)hLj-^L)/kBT _|_ I 

dE 



Similarly, the total energy is given by 



(Bl) 



+ 00 



(B5) 



E^.l- 



-f OO 



EdE 



The total number of bosons is given by 



1 



^ p((j + l/2)riu;-M)/fcBT _ 1 



1 



^^{huj /2 — fi) /Re 



-\-oo 



^ f,{(3 + l/2)nuj-iJi.)/kBT _ I 



(B2) 



(B3) 



The second summation term can be evaluated by integral 
using the Euler-Maclaurin formula: 



Running two independent indexes is computationally 
very inefficient and slow. We therefore assume loi ~ uj 
and UJ2 = kuj where A: is a positive integer. The summa- 
tion is run until (i + kj) = Q which, with Q large enough, 
yields a good approximation. Now N"^^ becomes: 



i+k]=Q 
N2,= E 

Q 

= E 



1 



4{i+k])tu^ + {l/2+k/2)hui-fj.]/kBT _|_ I 



floor{i/k) + 1 



(B6) 



i=0 



'irv^ + {l/2+k/2)hui-ti\/kBT _|_ ^ ' 



b 



E(|y![^^''"'(^)-^'''"H«)] + i? 



where floor{j/k) is the nearest integer less than or equal 
to j/fc. In this way only one index is present, allowing to 
improve significantly the computational speed. 



where B2 = 1/6, = -1/30, ...are the Bernoulli num- 
bers, and R is the remainder term. In the case of 
/cbT ^ fiw, the first term (fc = 1) in the expansion is 
sufficient, giving: 
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p{nuj/2-fj,}/kBT _ I 
+°° dx 



^ ^[(x + l/2)h^-ti\/kBT _ I 
1 

2 \(,{3huj/2-ti.)/kBT „ J^l 
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12KbT ^g(3fia;/2-/x)/fcBT _ ly 
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